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We present an analysis of the Korringa ratio in a dirty metal, emphasizing the case where a 
Stoner enhancement of the uniform susceptibilty is present. We find that the relaxation rates are 
significantly enhanced by disorder, and that the inverse problem of determining the bare density of 
states from a study of the change of the Knight shift and relaxation rates with some parameter, 
such as pressure, has rather constrained solutions, with the disorder playing an important role. 
, Some preliminary applications to the case of chemical substitution in the Rb3_a;Ki,C6o family of 

■ superconductors is presented and some other relevant systems are mentioned. 

INTRODUCTION 

Nuclear magnetic resonance (NMR) is an important 
. experimental tool for elucidating the physics of electronic 
materials whose properties are determined by strong 
, electron-electron correlations. For example, NMR data 
for high- Tc superconductors and for low dimensional con- 
, ductors has given us important information about their 
• electronic properties. Furthermore, such data are key 
I benchmarks against which theoretical approaches must 
, be measured. 

■ Basic NMR measurements in a metal involve the mag- 
, netic resonance field shift K and the spin relaxation rate 

d • 1/Ti of a given nucleus in the material. Q The effects 
^ ] of electron-electron interactions on these quantities have 
I ' been studied theoretically for nearly forty years [^-^ and 
' O I it is often taken for granted that a deviation of the "Ko- 
^ rringa ratio" /C = l/iTiTK"^) from a constant is a sig- 
^ nature of the effect of strong correlations. In particular, 
• • , the standard canon is that if the Korringa ratio exceeds 
unity, it is evidence for an antiferromagnetic enhance- 
ment of the electron susceptibility away from zero wave 
5^] ■ vector which enters 1 /Ti , |^ while if IC is less than unity it 
is because there is a (ferromagnetic) Stoner enhancement 
of the uniform susceptibility which enters K. |^,^ It is 
our purpose to point out that the well-known corrections 
to the susceptibility due to disorder [|j influence the Kor- 
ringa ratio in an important way which must be considered 
when drawing conclusions from the experimental results. 
We shall show that disorder significantly enhances the 
Korringa ratio and that it must be included when an 
attempt is made to estimate Stoner factors from IC. In 
what follows, we shall not consider the case of antifer- 
romagnetic enhancements of the wave- vector dependent 
electron susceptibility. Our discussion of the ferromag- 
netic Stoner enhancement is sufficient to illustrate the 
point. 

When the resonance field shift is determined solely by 
the local spin polarization of the conduction electrons, it 
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is called the Knight shift and it is proportional to the hy- 
perfine coupling gjN^ and the uniform conduction elec- 
tron susceptibility x(q ^ 0, w = 0). Here, 7Ar is the 
nuclear gyromagnetic ratio so that A has dimension of 
magnetic field. At the same time, 1/Ti depends on 
and on a density of electron magnetic excitations which is 
given by the imaginary part of x- The latter is averaged 
over all wave vectors weighted by a form factor ^ which 
depends on the position of the nucleus in the unit cell 
and may also include the effect of a transferred hyperfine 
interaction. Thus we have 

K = AxiO,0)/^iB (1) 

l/Ti - kBT{j/fiB)' ^'(q)x"(q: ^o)/cJo, (2) 
q 

where A{q) is the "q-dependent hyperfine" coupling 
which involves the form factor mentioned above. 

The case of free electrons is particularly simple. The 
susceptibility is of the Lindhard form: 

x'lIq, cj) = Xp[1 - iq/2kF f/3], xliq, ^) = {TTXp/2){u;/vFq) e{vFq - uj), 

(3) 

to leading order in q,uj. Here xp = 'rn*kpfi'^/7r^ = 
2/3oMb is the Pauli susceptibility, po the one-spin den- 
sity of states at the Fermi level. With one atom per unit 
cell, the form factor is unity and Eqs. (1, 2) combine to 
give the well-known Korringa relation: 

l/TiT = ^rkBhN/^iBfKyA (4) 

which makes it appear that 1/Ti is proportional to the 
square of the uniform susceptibility, or to the square of 
the density of states. Conclusions have often been drawn 
from this. However, this dependence is accidental; from 
Eq. (2) we see that l/Ti oc x/r, where F is a charac- 
teristic energy associated with the susceptibility at the 
important range of wave vectors in the sum. In any case, 
the Korringa product will cease to be temperature inde- 
pendent whenever the susceptibility is temperature de- 
pendent. 

In what follows, we show that diffusive corrections to 
x(q, uj) which occur in a disordered metal can also 
enhance l/Ti and must be considered when the mean 
free path is short, as is the case, for example, in some 
high- Tc samples and certainly in MaCeo- |@1 

THE CALCULATION 

We start with the usual RPA result for the suscepti- 
bility x(g, uj) 

xiq^Lu)-' ^Xo{q,co)-'-U/2, (5) 

where U is the effective Coulomb interaction and xo is 
the dynamical susceptibility, incorporating the effects of 
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disorder, but not those of interactions. We take units 
such that h = fiB = ^- We can express xo as a sum over 
impurity Hues, renormahzing propagators and vertices to 
leading order in the impurity scattering rate. This gives 
rise to the usual diffusion propagator |^ for small q,(jj. 
We can express the result in the form 

Xdii,^) -XL{q,uj)—^ — — , (6) 

where xl is the Lindhard susceptibility for the undisor- 
dered, non-interacting system. Dq = ^vpX is the dif- 
fusion constant with A the mean free path. To lowest 
order in disorder and interaction, the mean free path is 
the same one as enters the resistivity (for point impuri- 
ties). However, in general, the spin-triplet electron-hole 
channel mean free path enters the susceptibility, and it 
can be quite different (smaller) from the one involved in 
charge transport in the interacting system We have 
constructed Eq. (6) so that it has the correct behavior 
for both large and small q. Thus, it is correct to lead- 
ing order in the disorder, in the limit of low frequencies, 
and to O(g^) in the q dependence. We note the leading 
behaviour for small to 

x'^{q,uj)/uj = po{-^ + -^), (7) 
2qvF -Do? 

This expression is valid for essentially all the range q < 
2kF, and shows that the Landau damping [the first term 
in Eq. (^] ceases to dominate when the dimensionless 
number TTqX/6 < 1. Then the diffusion mechanism domi- 
nates the damping. Also to lowest order, the real part of 
X is unchanged by disorder and is still given by Eq. (^. 

In order to calculate the NMR relaxation rate, we need 
x"('?:Wo)/wo which is obtained from Eq. (^) as 

x"{q,^o) ^ Xo(g,^o) 1 

LOO [l-C/x;,(<Z,0)/2]2- ^> 

In Eq. (^), we neglected a term [f/xp (g, wq)]^ in the de- 
nominator since at the low NMR frequency cjqj it van- 
ishes as Wg. The Stoner denominator in Eq. (8) is sen- 
sitive to Xo(9jO)- Instead of using the approximation of 
Eq. (6) for it, we can find the result which is valid for all 
q as follows: 

At w = there are no diffuson (i.e. vertex) corrections 
to Xo(95^) and it is given simply by 

Xo(q,0) = -2 5:!tol^. (9) 

When the impurity scattering rate F — vf/2X is larger 
than the temperature (we take h ^ 1), 

n(k) = arctan (10) 

where the kinetic energy ek is measured from the chemi- 
cal potential. 
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We perform the angular integrals in Eq. (g) in two and 
three dimensions and find 

(2D) x^(,,0) = 2^^^dx-^ (11) 

/>oo 

(3D) XoiliO) ~ — / xdxn{xkp)\n 
^ Jo 



t + x 



t — X 



(12) 



where x = k/kp and t — q/2kp^ rj is kp\, and po is the 
one-spin Fermi surface density of states in the appropri- 
ate dimension. These may be evaluated in closed form as 
g ^ 0, where we find x'(0,0) = 2pon(k = 0) in 2D and 

X'(0,0) = 2po\/l + y(l + ?7-2)/y2 in 3D. These imply 
that disorder changes the effective density of states at the 
Fermi level by a small amount given by these formulae, 
and for consistency, we must use these factors in esti- 
mating the enhancement of the susceptibility or Knight 
shift. 

The result for the NMR relaxation rate is obtained by 
using Eqs. (7, 11, 12) in Eq. (8) and integrating over q. 
The results are conveniently expressed in terms of the 
free electron result: 

l/TiT = TTkB{iNPoA/2f ■ S{pQU,kpX), (13) 

where the enhancement factor S{poU, kpX) contains all 
the effects of interactions and disorder. In terms of the 
reduced variables, we have 



(2D) S{poU,7j)= [ dt 

(3D) S{poU, r;) = 2 / tdt 
Jo 



1 + {2/Tn]t) 
[l-Ux'o{t,OW 

1 + {3/TTJ]t) 

^-^^[i-Ux'oitM' 



(14) 



(15) 



In Eq. (15), we have introduced the infra-red cutoff k 
which is determined by an inelastic scattering ("Thou- 
less") length beyond which the diffusion ceases. In par- 
ticular, K ~ 1/ y/rjkpXi^ where Xi is the inelastic mean 
free path. These integrals can be evaluated for U — Q, 
where we find S{poU,ri) = 1-1- 6/{'Krj) in 3D and in 2D, 
S{poU,ri) = 1-1- \og{kpXXi)/{TTr]). Note that S is in gen- 
eral bigger than unity due simply to the enhanced density 
of states at low energies implied by the diffusive character 
of spin fluctuations in a disordered metal. 

A leading approximation for S is obtained from Eqs. 
(14,15) by replacing Xo(^iO) 3D, the approximate 

form 2po(l — i^/3) of Eq. (3). In 2D, the same approx- 
imation is XQ(t,0) ~ 2pQ. Therefore, the approximate 
expressions for the enhancement factors are 

(2D) S{poU,v)= I (16) 

(1 - Upo)^ ITT] 



(3D) 5(pof/,?y) ^ ^ 



[l-Upof U + y^ TTkpX 



1 + y V 



1 

arctan y 



, (17) 



where y = a/J7po/[3(1 - U po)]- 
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We can also calculate the dimensionless Korringa ratio 
K. (essentially the ratio of l/TiT to the square of the 
Knight shift K ) as 

lC = S{x^lxf- (18) 

This number is also a measure of the interactions and dis- 
order present in a Fermi system, and is of course merely 
an alternative description to that implied by S. How- 
ever, in some systems such as M3C60, the Knight shift 
is not easy to estimate, since the measured shift has to 
be apportioned into the chemical and Knight shifts with 
large uncertainties in both, whence we prefer to present 
both K, and S. 



CONCLUSIONS 

We show some figures which illustrate the discussion 
of this paper. In Fig. 1, we show how disorder affects the 
enhancement factor S in three dimensions. Even when 
the interaction (Stoner) enhancement is small, that is 
C/po < 1, iS can be enhanced by a factor 10 when the 
disorder is sufficiently great {kpX — 1). Fig. 2 shows 
the same situation in two-dimensions. The insets in the 
above two figures show that the approximate formulas of 
Eqs. (17, 18) for S in three and two dimensions are suffi- 
ciently accurate except at the largest values of the Stoner 
factor. Since experiments usually give the Korringa ratio 
K, directly, we plot it in Figs. 3, 4 as a function of disor- 
der and interaction strength. In all the two-dimensional 
plots we chose the cutoff k = 1/ -ySM^A 

In Fig. 5, we show a result which is potentially useful 
for a controlled study in MaCgo where properties appear 
to depend rather universally on lattice spacing. For ex- 
ample, in the series M = K3_j;Rbj, the variations in Tc, 
susceptibility |^ and l/Ti |^ have been measured. At 
present, there seems to be no substantial variation of re- 



sistivity with X, [g^ 11| thus we take kpX fixed in the 
plot which shows the change in 1/TiT relative to that of 
the uniform susceptibility x ^ the bare density of states 
Pq is varied. At the present stage, the uncertainties in 
the measurements prevent us from making a good infer- 
ence of the value of the Stoner enhancement or of the 
bare density of states, but we expect that with increased 
experimental accuracy, we should be able to exploit the 
present method to infer such quantities, which are of fun- 
damental interest. 

As a matter of fact, Knight shift and NMR relaxation 
rate data for dilute non-magnetic alloys have been used 
for man y y ears to study electronic structure in these ma- 
terials. [iSfl In view of the disorder enhancement to l/Ti 
which we have discussed, the systematics of the behav- 
ior with impurity concentration should be included in all 
such analyses. More recently, NMR has been used to 
study the properties of the rare-earth and actinide heavy 
electron materials, [ p^ For example, the superconducting 
compounds Ui_a;Tha;Bei3 are of special interest because 
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of the dramatic effects caused by the Th doping. There 
is a large increase in 1/Ti just above Tc as x goes from 
to 0.033 which is not well-understood but has been 
attributed to a density of states effect. Again, it is impor- 
tant to include the effects of disorder on the relaxation 
rate, as we have discussed. 

A system that does show the su gge sted enhancements 
in l/TiT, is the well-studied Si:P where the ob- 

served NMR relaxation rates are approximately three or- 
ders of magnitude larger than the clean free electron val- 
ues. Here there is the added complication that there are 
two species of relaxing fluids, the dirty electron gas as in 
the present paper, and the singlet pairs formed by local 
moments that fluctuate and relax the nuclei in a distinc- 
tive fashion as described in [0. The detailed quantita- 
tive breakup between these terms is not easy to do and is 
not attempted here. The point we wish to make is that 
the free electron gas component contributes a very large 
amount (see Figs, f-2) to the observed 1/TiT when both 
Stoner enhancement and disorder are present. 
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FIG. 1. S versus X = poU in 3-dimensions for 
kpX = 1, 3, 5, 10, oo from top to bottom. Inset compares the 
leading approximation (dashed lines) Eq. (18) with the exact 
numerics (solid lines) for A;fA = 1, 10 

FIG. 2. S versus X = poU in 2-dimensions for 
fcfA = 1, 3, oo from top to bottom. Inset compares the lead- 
ing approximation (dashed lines) Eq. (17) with the exact 
numerics (solid lines) for /cfA = 1, 10 

FIG. 3. The Korringa ratio JC versus X = poU inS- dimen- 
sions for fcfA = 1, 3, 5, 10, oo from top to bottom. 

FIG. 4. The Korringa ratio /C versus X = poll in 2- dimen- 
sions for ^fA = 1, 3, 5, 10, oo from top to bottom. 

FIG. 5. The logarithmic 

derivative P = d\n[{poUfS{p»U,kF\)]/d\n[p»U/{l - poU)] 
versus X = poll in 3-D. This corresponds to a variation of the 
bare density of states for a fixed U and the resulting change in 
l/TiT relative to that in the uniform susceptibilty x on a log- 
arithmic scale, i.e. dln(l/rir)/dlnx. Upper (lower) curve is 
for krX = l(oo) 
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